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ABSTRACT 

The  Paarbicle-in-Cell  method  for  multidimensional  fluid- dynamics  cal¬ 
culations  has  been  used  on  the  IBM  Stretch  Computer  to  study  the  motion  of 

a  hot  gas  bubble  created  by  a  nuclear  explosion  in  the  earth's  atmosphere, 

20 

In  the  example,  a  total  energy  of  8.24  x  10  ergs  was  deposited  at  an  al¬ 
titude  of  95  km.  After  2.85  sec  the  originally-heated  air  was  moving  up¬ 
wards  with  a  mass-avereiged  velocity  of  5*5  km/ sec  and  was  still  accelera¬ 
ting  somewhat,  probably  to  an  asymptotic  value  of  4,0  km/sec.  Along  the 
cylindrical  axis  at  2,85  sec,  the  upwards  velocities  varied  from  0.0  km/ sec 
at  the  bottom  to  5,8  km/sec  at  the  top;  most  of  the  central  material  was 
moving  at  velocities  greater  than  4,0  km/ sec.  This  report  also  includes 
details  concerning  bubble  shape,  energy  histories,  and  shock  configurations. 
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INTRODUCTION 

The  deposit  of  considerable  energy  into  a  i^herical  section  of  the 
earth's  atmosphere  would  be  followed  by  subsequent  non-spherical  air  move¬ 
ments,  with  the  initially  heated  region  elongating  vertically  and  rising. 
We  have  studied  in  detail  a  particular  example  and  present  the  results  be¬ 
cause  of  several  features  which  apply  generally  to  this  class  of  phenomena. 
The  study  was  purely  theoretical;  the  fluid-dynamics  equations  were  solved 
numerically  by  a  hi^- speed  computer,  and  no  direct  comparisons  with  ex¬ 
periments  are  presented. 


THE  SPECIFIC  PROBLEM 


The  model  is  a  simplified  one  of  the  conditions  which  mi^t  be  fovind 
shortly  after  a  hi^-altitude  nuclear  explosion.  The  specific  example  was 
suggested  by  C.  Longmire,  At  initial  time,  t  =  0,  the  disturbed  air  was 
a  heated  sphere  with  center  at  95  Ion  above  ground  and  radius  of  2,7  km. 

It  was  assumed  that  no  motion  had  occurred  but  that  the  ^here  of  air  had 
been  heated  to  a  uniform  specific  internal  energy  of  800  (km/ sec)  and 
had  a  radial  velocity  whose  magnitude  was  given  by 
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Here  Z  is  the  hei^t  above  ground,  vdth  being  the  center  height.  The 
atmospheric  density  variation  was  taken  to  be 


p=  1.2x10  exp 


z  -  V 

6  km  _ 


Finally,  the  equation  of  state  of  the  air  was  taken  as  polytropic,  with 
the  ^ecific  heat  ratio  in  the  heated  ^here  being  y  =  1.6667,  while  that 
outside  was  7  =  1.5555*  These  values  were  chosen  so  as  to  account  for 
dissociation  and  ionization.  The  acceleration  of  gravity  was  considered 
negligible  except  insofar  as  it  created  the  initial  density  distribution. 


THE  METHOD  OF  SOEUTION 

Since  the  problem  involves  two  space  dimensions  (R  and  Z  in  cylin¬ 
drical  coordinates)  and  results  in  considerable  distortion  of  the  air,  we 
have  used  the  Particle-in-Cell  (PIC)  method  as  currently  written  for  the 

IBM  Stretch  Computer  in  a  Group  T-5  code  called  DUX.  The  basic  technique 

* 

has  been  described  in  detail  in  IA-2501  and  the  references  mentioned 
there.  DUX,  however,  uses  several  sU^t  modifications,  particularly  of 
importance  in  the  treatment  of  the  energy  equation.  Details  are  avail¬ 
able  from  the  authors;  thus  we  give  only  a  brief  description  here. 

The  computing  process  employs  basically  an  Eulerian  viewpoint. 

Thus,  a  cylindrical  region  of  ^ace  about  the  originally- heated  inhere 
was  chosen  to  be  large  enough  to  contain  the  region  of  interest 

*F.  Harlow,  IA-250I,  April  1959. 
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througiiout  all  times  of  significance.  In  this  case,  the  cylindrical  re¬ 
gion  extended  from  12  km  below  the  center  to  58,4  km  above,  and  had  a 
radius  of  19,2  km.  A  slice  through  the  cylindrical  region  was  then  sub¬ 
divided  into  a  finite-difference  mesh  of  cells  whose  sides  were  0,6  km 
long.  The  air  was  represented  by  particles  in  the  cells,  whose  motions 
gave  the  configuration  changes.  One-hundred  per  cell  were  used  in  the 
originally-heated  air,  while  four  per  cell  were  outside, 

A  PIC-method  calculation  proceeds  through  a  sequence  of  finite  time 
steps  in  each  of  which  the  fields  of  velocity,  energy,  position,  etc.  are 
changed  by  sm^n  amounts.  The  results  are  observed  through  summary 
prints  given  by  the  computer  and  by  configuration  pictirres  given  through 
the  SC-4020  Microfilm  Recorder.  From  the  summary  prints  there  can  be  ob¬ 
tained  the  time  histories  of  such  functionals  as  the  totals  of  momentima, 
kinetic  and  internal  energies,  and  volumes.  Detailed  profiles  of  the 
local  densities  of  mass,  energy,  etc,  can  also  be  obtained  from  the  com¬ 
puter  whenever  desired.  From  a.1 1  of  the  available  data,  the  conclusions 
of  significance  can  then  be  drawn, 

RESUITS  OF  THE  CALCULATION 

The  configurations  are  shown  in  Fig,  1  a-g,  in  which  the  unretouched 
dot  positions  are  shown  for  a  sequence  of  times.  In  the  last  picture,  for 
which  the  time  is  t  =  2.85  sec,  it  can  be  seen  that  the  shock  pulse  has 
broken  away  from  the  "piston"  of  initially-heated  air,  and  is  followed  by 
a  rarefaction.  Vertical  stratification  of  the  air  is  caused  partly  by 
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shock-rarefaction  interactions,  but  may  also  arise  because  of  the 
initially-coarse  resolution  by  the  finite-difference  mesh  (whose  cells 
are  not  shown). 

Energy  histories  are  shown  in  Fig.  2.  Noteworthy  is  the  nearly- 

constant  behavior  of  external  kinetic  energy  (Curve  C).  The  vertical 

momentum  of  the  initially-heated  air  is  shown  as  a  function  of  time  in 

/  8  \ 

Fig.  5.  Division  of  this  by  the  total  mass  in  the  bubble  (1,02  X  10  gm) 
gives  the  mass-average  velocity  which  can  be  read  from  the  same  figure. 

Finally,  Fig.  h  shows  a  profile  of  the  vertical  velocity  along  the 
cylindrical  axis  at  a  time  t  =  2,85  sec.  The  dat\mi  points  are  cell- 
wise  values  (which  have  the  PIC-method  fluctuations  to  be  escpected  for 
this  situation),  while  the  line  is  an  eye-drawn  mean. 

While  this  exainple  refers  to  a  specific  initial  configuration,  the 
qualitative  natxnre  of  the  results  would  apply  to  a  variety  of  similar 
problems.  PIC-method  calculations  would  give  details  in  cases  where 
quantitative  details  are  required. 
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Figure  la.  Sequence  of  particle  configurations  shoving  the  air  notions 

about  the  initially-heated  bubble.  The  squares  are  four  kilo¬ 
meters  on  a  side,  t  =  0.0  sec. 
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Figure  Tb.  Sequence  of  particle  configurations  showing  the  air  motions 

about  the  initially-heated  bubble.  The  squares  are  four  kilo¬ 
meters  on  a  side,  t  =  0.45  sec. 
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Figure  1c.  Sequence  of  particle  configurations  showing  the  air  motions 

about  the  initially-heated  bubble.  The  squares  are  four  kilo 
meters  on  a  side,  t  =  O.87  sec. 
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Figure  Id.  Sequence  of  particle  configurations  showing  the  air  motions 

about  the  initially-heated  bubble.  The  squares  are  four  kilo¬ 
meters  on  a  side,  t  =  1.3^  sec. 
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Figure  1e.  Sequence  of  particle  configurations  showing  the  air  motions 

about  the  initially-heated  bubble.  The  squares  are  four  kilo¬ 
meters  on  a  side,  t  =  1.75  sec. 


-n- 


Figure  If.  Sequence  of  particle  configurations  showing  the  air  motions 

about  the  initially-heated  bubble.  The  squares  are  four  kilo¬ 
meters  on  a  side,  t  =  2.41  sec. 
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Figure  1g.  Sequence  of  particle  configurations  showing  the  siir  motions 

about  the  initially-heated  bubble.  The  squares  are  four  kilo¬ 
meters  on  a  side,  t  =  2.85  sec. 
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Figure  2,  Energy  histories  for  the  bubble  and  surrounding  air.  A  and  B 
refer  to  the  initially-heated  air  and  are,  respectively,  the 
kinetic  energy  (xIO"''®)  and  internal  energy  (xlO'^y).  c  and 
D  are,  respectively,  the  kinetic  energy  (xlO"'^)  and  internal 
energy  (xlO-^^)  of  the  not-initially-heated  air. 
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Figure  4 


Velocity  profile  along  cylindrical  axis  for  a  time  t  =  2 
sec.  Dots  show  cell  values  while  line  is  eye-drawn  mean 


